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Abstract. Large-amplitude Alfven waves are ubiquitous in space plasmas and a 
main component of magnetohydrodynamic (MHD) turbulence in the heliosphere. 
As pump waves, they are prone to parametric instability by which they can generate 
cyclotron and acoustic daughter waves. Here, we revisit a related process within 
the framework of the multi-fluid equations for a plasma consisting of many species. 
The nonlinear coupling of the Alfven wave to acoustic waves is studied, and a 
set of compressive and coupled wave equations for the transverse magnetic field 
and longitudinal electric field is derived for waves propagating along the mean-field 
direction. It turns out that slightly compressive Alfven waves exert, through induced 
gyro-radius and kinetic-energy modulations, an electromotive force on the particles 
in association with a longitudinal electric field, which has a potential that is given 
by the gradient of the transverse kinetic energy of the particles gyrating about the 
mean field. This in turn drives electric fiuctuations (sound and ion-acoustic waves) 
along the mean magnetic field, which can nonlinearly react back on the transverse 
magnetic field. Mutually coupled Alfven-cyclotron-acoustic waves are thus excited, 
a nonlinear process that can drive a cascade of wave energy in the plasma, and may 
generate compressive microturbulencc. These driven electric fiuctuations might have 
consequences for the dissipation of MHD turbulence and, thus, for the heating and 
acceleration of particles in the solar wind. 



1. Introduction 

Large-amplitude Alfven waves are ubiquitous in space plasmas, and particularly 
prominent in the solar wind (Tu and Marsch, 1995; Bruno and Carbone, 2005). 
They are an essential component of magnetohydrodynamic (MHD) turbulence in 
the heliosphere and known to originate mainly in the solar coronal holes (Cran- 
mer, 2009). As has been shown in the ample literature, an Alfven mother (pump) 
wave is prone to parametric instability (Stenfio, 1976; Derby, 1978; Goldstein, 1978; 
Longtin and Sonnerup, 1986; Brodin and Stenfio, 1988; Hollweg, 1994; Wong and 
Goldstein, 1986; Vinas and Goldstein, 1991a, b; Stenfio and Shukla, 2000; Ruderman 
and Simpson, 2004; Stenfio and Shukla, 2007) by which it can generate cyclotron 
and acoustic daughter waves that may undergo kinetic effects (Araneda, 1998) and 
collisionless Landau damping (Inhester, 1990; Araneda et al., 2007). The continu- 
ous and wide interest in these waves also comes from their astounding properties, 
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namely that Alfven-cyclotron waves, like parallel magnetosonic-whistler waves, are 
nonlinear eigenmodes (Sonnerup and Su, 1967; Stenflo, 1976) of the MHD, and 
multi-fluid equations as shown below, for propagation along the mean magnetic 
field. 

Nonlinearly excited (Spangler, 1989) acoustic waves appear to be common in 
space plasmas as well, and density fluctuations (Tu and Marsch, 1995; Bruno and 
Carbone, 2005) are observed everywhere in the solar wind, although at a compar- 
atively low fluctuation level of merely a few percent. However, since compressive 
fluctuations can be clamped through kinetic effects, like Landau damping on the 
thermal ions and electrons, they can provide an effective dissipation mechanism 
for the nonlinear damping (Medvedev et al., 1997) of Alfven waves. Consequently, 
the understanding of the coupling between Alfvcnic wave activity and density or 
charge-density fluctuations is of paramount interest and importance in basic plasma 
physics, but alike in its applications to nonlinear processes in space (Stenflo and 
Shukla, 2007) and astrophysical plasmas. 

As will be shown in this paper, a coupled set of nonlinear second-order wave 
equations for the transverse magnetic fleld, the transverse gyromotion of any parti- 
cle species in the multicomponent plasma considered, and the related longitudinal 
electric fleld can be derived, which together describe the wave-wave interactions 
and their mutual forcing. These equations provide a physically and intuitively clear 
picture of the field and particle/plasma dynamics and allow us to understand the 
results of recent hybrid simulations of the parametric decay of Alfven waves and 
their effects on the plasma particles better. The main aim of this work is to provide 
algebraic derivations and physical explanations. A numerical trciatment of the full 
equations to be derived subsequently appears promising, yet is beyond the scope of 
this work. 

In analytical (Araneda et al., 2007), hybrid-simulation, and other numerical sim- 
ulation (Araneda et al., 2008) studies of the parametric instabilities of Alfven- 
cyclotron waves, it became obvious that ion trapping (Araneda et al., 2008, 2009) 
in the nonlinearly driven ion-acoustic waves and pitch-angle scattering by the trans- 
verse daughter waves were found to cause anisotropic heating of the proton core 
velocity distribution, and simultaneously to create a proton beam along the mean 
field (Araneda et al., 2008; Valentini and Veltri, 2009). These numerical results are 
in close agreement with observed kinetic features in the solar wind and support 
the observation that pitch-angle scattering (Heuer and Marsch, 2007; Marsch and 
Tu, 2001) is the key to understand the kinetic characteristics of thermal solar-wind 
protons. But only recently convincing evidence has been found for ion-cyclotron 
waves (Jian et al., 2009) to exist in the solar wind. Also simulations of electric field 
spectra (Valentini et al., 2008) have shown that the short-scale termination of so- 
lar wind turbulence is characterized by the occurrence of longitudinal electrostatic 
fluctuations. The spectra thus obtained seem to be consistent with the electrostatic 
waves actually measured in the solar wind (Bale et al., 2005) close to the Earth's 
bow shock, and in particular in the ion-cyclotron range (Kellogg et al., 2006). 

The present study will provide the foundation for insight into and further study of 
the processes occurring at macroscopic and microscopic scales in solar-wind turbu- 
lence, and thus will throw light on the related dissipation processes through kinetic 
cascades and wave-particle interactions (Marsch, 2006). The nonlinear equations 
derived here are used to describe an elliptically polarized Alfven wave as a simple 
but nontrivial example of their application. 
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— • (njV,) = 0, (2.1) 



2. The multi-fluid equations in conservation form 

2.1. Fluid equations in the wave frame 

In this section, we shall first recapitulate the basic multi-fluid equations (Goossens, 
2003) for a plasma consisting of electrons and various ionic species. We start from 
the fundamental conservation laws, and then try to make no approximations with 
respect to the field amplitudes in order to be able to discuss and analyze nonlinear 
waves and convected wave-like structures. We are going to use coordinates in the 
frame of reference moving with the wave, which has a normal to its front denoted 
by n and a propagation speed V = Vn in the inertial frame, or center of momen- 
tum frame that is defined below. This unit vector obeys the relation = 1. The 
coordinate in this moving frame is ^ = x — Vt, and all variables are assumed to de- 
pend on space and time only through ^, and thus spatial and temporal derivatives 
in the wave frame are reduced to derivatives with respect to ^. Such a coordi- 
nate transformation has been used by many authors, for example to study solitary 
waves in multi-ion plasmas (Hackenberg et al., 1998), or their stability properties 
(McKenzie et al., 1993). Therefore, by using comoving coordinates. Maxwell's par- 
tial differential equations in space and time, and similarly the fluid equations for 
the different species, can be reduced to simpler differential equations in terms of ^. 
The continuity equations thus reads 

d_ 

with Vj being the flow velocity of species j in the moving frame Vj = Uj — V, and 
rij is its number density. With the charge denoted as qj, we get the total charge 
density as 

j j 

which must obey Gauf5' law: 

47^(7 = ^ • E. (2.3) 

Similarly, the total current density is given by 

J = Y.qjnjUj, (2.4) 
j 

which has to obey Ampere's law: 

The second term of Eq. (2.5) is the displacement current. The conduction minus 
convection current density may be written 

j = J-aV = J2<^jVj. (2.6) 

i 

To complete the set of Maxwell's equations, we quote that the magnetic field must 
be divergence free 

^•S' = 0, (2.7) 
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and that Faraday's induction equation requires the curl of the electric field in the 
wave frame to vanish: 

^xE' = 0, (2.8) 

with the primed variables being defined in the wave frame. The Lorentz transfor- 
mation has been used to derive Eq. (2.8), since it gives the connection between the 

electromagnetic fields (c is the speed of light in vacuo) in the plasma's center of 
momentum frame and the moving wave frame through the relation 

E' = E+-V X B. (2.9) 
c 

Of course, the magnetic field remains invariant to lowest order in V/c, and thus 
B' = B. For later purposes, we define the mass density of particle kind j as 
Qj = rijUij, with the total mass density 

Q = Y.Qr (2.10) 

3 

The center of momentum velocity (for which we are free to choose C/ = 0) is defined 
as follows: 

qU = Y^q,U,. (2.11) 

Since we are interested in the individual ion and electron dynamics, we do not sum 
their momentum equations like in MHD, but use instead the separate multi-fluid 
equations. The individual momentum equation of species j can conveniently be 
quoted in conservation form in the moving frame, reading 

d_ 

The expression VjVj means a tensor in dyadic notation, and 1 the unit dyade. 

For the equation of the partial pressure, we may take a simple polytropic equation 
of state for the purpose of closure, and thus write 

Pj=Pjo(^Y\ (2.13) 

with some constant reference density, njo, and pressure, pjo. Equivalently, we may 
consider the (polytropic, with index 7j) entropy equation 

d_ 

The set of Eqs. (2.1), (2.3), (2.5), (2.7), (2.8), (2.9), (2.12), and (2.13) is closed and 
sufficient to calculate all the independent but coupled variables. In what follows we 
shall assume a reduced geometry. 

2.2. Reduced multi-fluid equations in one- dimensional geometry 

We consider a one-dimensional spatial setup, i.e. a dependence on only one spatial 
coordinate and take components with respect to n, where the unit vector may 
correspond to the wave unit vector fe in Fourier variables. Thus we generally define 
the components 

^ = n • (a; - Vt), (2.15) 



{mjnjVjVj +p,l) = qjnj (^E' + x B'^ . (2.12) 



^r-^ln{pj0r')=O- (2-14) 
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Vj ^Vjnfi + Vjt, (2.16) 
Vjt = il-nn)-V,. (2.17) 

Transverse components are obtained by projection perpendicular to the longitudinal 
direction. The corresponding magnetic field components are defined as = n ■ B 
and = (l — nn) • -B, from which it follows that n x Bt = fi x B, and that 
Bt-h = Q. The fluid equations then read as follows. For the longitudinal momentum 
conservation, we have 

i(„,r„V5. + ^) = S!i^(K + i(V.xB0.*.), (2.18, 

and for transverse momentum conservation we obtain the equation 

d , , qjrij /„, 1 



_ (n,y,„F,0 ='^[^E[ + - {V, X B')t j • (2-19) 
The longitudinal magnetic field is strictly conserved and thus constant: 

^=0, (2.20) 
and the charge density obeys Gaufi' law 

^=4nJ2q^n,. (2.21) 

3 

If we take after Eq. (2.8) the curl of E' in the wave frame, we find 

, d . d£;; . dEt . /I,,. dBt\ 

n— X E' = hx = = n X — ^ + n x -Vh x — ^ , 2.22 
dC d^ d^ \c d^ J 

where the last part of the equation follows from the previous Eq. (2.9). In conclusion, 
E[ is constant and can be set equal to zero. Using this result, the transverse electric 
field is obtained through Eq. (2.9), which is equivalent to writing 

Et = --V{h x Bt), (2.23) 

c 

and makes the transverse electric field a dependent auxiliary variable being fully 
determined by the transverse magnetic field. For the longitudinal electric field com- 
ponent in the wave frame, one has 

E^ = E^ + -h- {Vh x B), (2.24) 

c 

which yields E'^ = E^, which is to be used in Gaufi' law (2.21). For the one- 
dimcnsional spatial geometry chosen here, the longitudinal current density can be 
written as 

Since the curl of B has only transverse components, we obtain from Eq. (2.25) that 
the longitudinal current density in the wave frame must be strictly constant, since 
it is given by 
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where the individual particles fluxes {Fjn) are conserved according to the longitu- 
dinal continuity equation, which expresses flux conservation in the form 

^ in.Vj^) = ^ = 0- (2-27) 

The transverse component of Ampere's law including the induction current can be 
cast in the form 

whereby we note that Vjt = Ujt, since V has no transverse component, but Vj^ = 
Ujn — V. The right-hand side of Eq. (2.22) gives an expression for the gradient of 
the transverse electric field component. It can be inserted in Ampere's law, which 
thus can be written as 

If {y/cf <C 1, which we will assume in the remainder, then the displacement 
current term can be safely neglected. It is certainly needed if one wants to make the 
transition to free electromagnetic waves correctly, in which we are not interested 
here. We will therefore not keep this term anymore. The basic equation for the 
magnetic field, which is Ampere's law for the transverse component, then reads 

^ E = ^inxB^). (2.30) 

i 

The magnetic field is free of divergence, which means in our geometry and variables 
that Eq. (2.20) is fulfilled with B„ = B'^. The last two equations together fully 
determine the vector magnetic field, given the current density is provided. The 
transverse electric field is obtained from Eq. (2.23) and the longitudinal one from 
GauB' law Eq. (2.21) with = E'^. 

Next we quote again the transverse momentum equation for each species, which 
reads 

4 {njVj^Vj^) = -^nj {Vj xB\ = ^ {Vj^n xB, + V x fiB^) , (2.31) 

whereby the electric field term has been written out in detail. Similarly, the longi- 
tudinal momentum equation reads 

which must be supplemented to obtain closure by the entropy or pressure equation, 

±ln(p,g-''^)=0. (2.33) 

In what follows it turns out to be convenient to use the natural spatial and 
temporal scales of the multicomponent plasma, which depend on the various fluid 
and field parameters. We define a longitudinal gyration length for the species j as 

K„ K„ 1 , . , 
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which corresponds to the gyration radius calculated with the longitudinal velocity 
instead of the perpendicular one. It is implicitly dependent on nj via the drift speed 
and Eq. (2.27). Another interesting length is given by the strictly constant quantity 

' 4^9,^;„ '\V,^) r,\uj ^ ^ 

Here we have introduced the Alfvcn speed, V^j = B\l {ATmjnij), based on the mass 
density gj of species j only, and the respective plasma frequency uij = AnqjUj/nij 
and gyrofrequency flj = B^qj/{mjc), carrying the sign of the charge qj. Note 
that Lj is strictly constant, and its inverse sums over all species up to zero, i.e. 

1/Lj = 0, because of the condition of zero longitudinal total current: — 
aV = QjnjVjn = 0. The standard Alfven speed based on is obtained by the 
summation 

Concerning the compressive dynamics, it is important to note that longitudinal 

and transverse motions are coupled through pj and rj = Vj^/^j, i.e. through the 
particle number density, when the mass continuity equation (2.27) and entropy 
equation (2.33) are exploited. 



3. Wave equations 

In this section, we recast the basic equations for the fields and plasma multi-fluid 
parameters into the form of coupled wave equations. For that purpose, we shall not 
straightforwardly Fourier transform them but rather rewrite them, by use of mul- 
tiple differentiation, in such a form that we finally obtain single "wave equations" 
for the electric field and magnetic field components. Remember that in the moving 
frame all variables depend solely on the coordinate ^ = x fi — Vt. We first consider 
the pressure equation again. We may use for the species' sound speed the standard 
definition 

and can then re-evaluate the momentiim eqiiation by use of 

d f pj\ 2 



which together with Eq. (2.27) allows us to quote the longitudinal momentum 
equation in the form 

It is convenient to introduce the effective Debye length Xj of species j as follows: 

1 uj] 

the sum of which gives the total Debye length, still including the differential drifts: 

(3.5) 



-=y- 
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Both Xj and Ad are not necessarily real quantities. Each species brings in its own 
length scale Aj. It is also convenient to introduce the second-order wave operator 
(which parametrically depends still on V via the Vjn)- 

rP 1 

Finally, one obtains a driven wave equation for the longitudinal electric field 

-y- 



VeE^ = -Y.V2 i^jt xBt)-n, (3.7) 



in which the transverse particle motions and electromagnetic fields show up through 
a nonlinear electromotive force, which is the summed contribution of the longitu- 
dinal components of the Lorentz forces acting on each species. This driving force 
acting on En resembles a "convection" electric field. When being decoupled from 
the transverse plasma and field dynamics, the longitudinal electric field equation 
just describes free electrostatic oscillations, such as Langmuir and acoustic waves, 
as we will show later. 

Let us return to the transverse momentum equation and rewrite it by exploiting 
the mass continuity equation. Then it is straightforward to derive 

Using this equation, we can rewrite the normal component of the convection electric 
field which occurs in Eqs. (3.3) and (3.7) as follows: 

1 , _ , B„ ri d /' I r, \ rrii d / 1 o \ 

E,n . - (y,. X BO • n = -^f ^ (-n) - -^d^ (2^^0 ' ^'-'^ 

a relation which is going to be used later. If the module of the transverse plasma 
velocity of species j is constant, then its convection electric field Ejn vanishes. Using 
Ampere's law Eq. (2.30) and the previous Eq. (3.9), wc can derive by vector cross 
multiplication of Eq. (2.30) with n and subsequent scalar multiplication with Bt 
another conservation law: 




0jV%-JzBt\=O. (3.10) 



If the integration constant is zero, this equation expresses equipartition between 
the transverse total particle kinetic energy and the transverse magnetic energy, like 
it is the case in a classical MHD Alfven wave. 

We consider now the equation for the transverse magnetic field. By differentiation 
of Ampere's law, we can obtain a second-order nonlinear wave equation for the 
transverse magnetic field 

^nxB, = -^..(n.^ + V,.-^j. (3.11) 

Here, it is convenient to introduce the skin depth or inertial length ^j of species j 
defined as follows: 

1 -I 
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the sum of which gives the total skin depth: 

i=Ei' (3-13) 



& J 3 



where each species brings in its own length scale 1^ . It is again convenient to intro- 
duce a second-order wave operator 



d^ 1 



= 372 - 72 • (3-14) 



Using this, one finally obtains a driven wave equation for the transverse magnetic 
field 




(3.15) 



On the right-hand side, the transverse currents appear and the longitudinal charge 
density variations and related electrostatic effects show up through the nonlinear 
electromotive force, which involves the longitudinal electric field. Note that this 
nonlinear driver contains the natural length scales (density dependent) of all the 
species involved. Of course, this wave equation seems, without further approxima- 
tion, quite formal but elucidates the nature of the coupling of the unforced trans- 
verse magnetic field (dynamics described by the operator Db) with the compressive 
electrostatic fiuctuations and transverse plasma motions. When being decoupled 
from the plasma currents and electric-field (no charges) dynamics, this transverse 
magnetic field equation just describes the finite penetration of the magnetic field 
into the skin layer of the plasma and results in its exponential decline on the length 
scale ^s- To gain better insight into the terms contributing to Eq. (3.15), we may 
rewrite it also in the form 

T^bB, = -\^^ \ —Vj, + 2_ 2 (n X V,t) ] ■ (3.16) 



3 



Similarly, one can also rewrite the driven wave equation for the longitudinal electric 
field in the concise form 

VeE^ = Y,^Ej^, (3.17) 

3 ^ 

where we remind that Ejn can be derived after Eq. (3.9) from a potential that is 
given by the transverse kinetic energy of species j. These two coupled nonlinear 
equations are completed and closed by Eq. (3.3) for the density and Eq. (3.8) for 
the transverse velocity of each species. 

So far, we made neither any approximation nor linearization, but just inserted the 
original momentum equations (3.3) and (3.8) into the differentiated Ampere's and 
GauB' laws. Therefore, the above equations depend in a highly nonlinear manner 
still on the different number densities nj. Yet the transition in Eqs. (3.16) and (3.8) 
to the incompressible limit is simple, because then the electrostatic nonlinear forcing 
terms E^ and Ej^ vanish, and the plasma frequency coj and the inertial length ij, 
respectively, become constants as defined by the fixed background density of species 
j. We discuss this in the section after the next one. 
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4. Eigenmodes and driven waves 

4-1- Nonlinear Alfven-ion- cyclotron waves 

In order to maintain the linear form of the original equations, it is convenient to 
introduce new variables, relating to left- and right-hand polarized fields, which are 
defined as follows: 

B± =B,±{hx Bt) , (4.1) 

vt = Vt±(nx Vt) . (4.2) 

These field variables are orthogonal, i.e. Bf ■ Bf = 0, and • Vj^ = 0. By taking 
the cross product of Eqs. (2.30) and (3.8) with the unit vector n, we obtain after 
some algebra the equations of motion for the circular transverse variables: 

^-V'^Vl. (4.3) 



dC c ^ 

j 



dV^ 

jt 



de 'r,\B^-' -^v' ^^'^^ 

Let us now first consider nonlinear incompressible solutions. Since En = 0, conse- 
quently quasineutrality strictly holds, qjUj = 0. The velocity fields and magnetic 
field must have constant modules and be aligned, so that after Eq. (3.9) their re- 
spective vector cross product, and thus also Ejn vanishes. As Uj, rj, Vjn, cj, and 
ojj then all are constant, we can solve the resulting linear set Eqs. (4.3) and (4.4) 
by Fourier tansformation (FT), without putting any limitations on the amplitudes 
of Bf or V^, other than from Eq. (3.10) which implies that the magnetic field 
amplitude is also constant. As usually, FT means that d/d^ — > i fc, and thus we can 
invert the transverse momentum equation (3.8), which yields with the normalized 
wave vector, kj = kVjn/^j, the complex vector relation 

This result can be inserted into the FT of Eq. (3.16) to obtain the algebraic wave 
equation 

(k-^ + 1") Bt{k) = E ^13^ (^*^^^ + '^'^'^ ' (^-^^ 



which may also be written as 





^t(A:) = H V (n X B,{k)), (4.7) 



and which yields, by taking the vector cross product of Eq. (4.7) with n and by 
resolving the resulting two equations, the two dispersion relations describing left- 
and right-hand polarized waves as follows: 



with the fractional mass density gj = Qj/ Q. Eq. (4.8) is nothing else but the standard 
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dispersion relation for the Alfven/ion-cyclotron (and magnetosonic- whistler) waves 
in a multicomponent plasma with the differential drifts contained in Vjn and for 
parallel propagation (see e.g. Davidson, 1983), yet which applies here to arbitrarily 
large wave amplitudes. Apparently, the wave frequency is obtained by the Doppler 
shift formula, oj = kV, where V is hidden in Vjn = Ujn — V, i.e. in Kj = kVjn/^j- 
Once the wave vector k = k{V) is known, the wave frequency is obtained as a 
function of the phase speed V. 

Let us consider the long- wavelength limit of Eq. (4.8), which when being ex- 
panded to second order in Kj reads 

l = E^^(|f)'(±«.)(l±«.)' (4-9) 

where the first term of the sum vanishes, since (|) ^.n = because of the 
quasi-neutrality condition, i.e. cr = in Eq. (2.2), and the zero-longitudinal-current 
constraint (2.26). The second term then yields 

. {U^^-2UjnV + V') 
1 = 2^ y2 • (4-10) 

3 ^ 

As a vanishing bulk speed, U ~ 0, may bo assumed, the center-of-momentum 
condition means that ffjUjn = 0, and thus we can solve for the phase speed in 
the center of momentum frame and obtain 



V = ±Va 



\ 



U- 



This is the phase speed of an Alfven wave in a multicomponent plasma, including 
field-aligned drift motions, leading to a slowing down of the phase speed. 

We shall now derive the incompressible Alfven-cyclotron wave without resort to 
the FT, but recourse instead on Eq. (3.9). Since Ej^ = 0, each velocity vector and 
magnetic field must be aligned, which generally implies that Vjt = ajB^- This can 
be inserted in Eq. (3.8) to obtain 

dv,, _i fv,n ^ (4_^2) 



d^ Tj \Bnaj / 
Twofold differentiation yields the simple harmonic oscillator equation for the gyro- 



motion: 

^''^ Ot=0, (4.13) 



with the squared wave vector defined as 

Since all species must spatially oscillate in the same way, the wave vector must not 
depend on the index j, i.e. we can put kj = ±k, which yields with Kj = krj two 
possible solutions for the desired proportionality coefficient: 
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Knowing the coefficient aj, we can use it in the wave equation (3.16) without electric 
fields and get another harmonic oscillator equation: 

(- 



[—+q^]B, = 0, (4.16) 



where the squared wave vector q is an abbreviation for exactly the same sum as 
appearing on the right hand side of Eq. (4.8). Since all velocities and the magnetic 
field are aligned, the wave vector q must be equal to k, and thus we again obtain 
the same dispersion relation as in the previous Fourier analysis. Finally, the polar- 
ization relation (with the plus sign for incompressible Alfven-cyclotron and minus 
for magnetosonic-whistler waves) reads 

With this result, we can further evaluate the conservation law (3.10), and after 
some algebra obtain the result that 

E^.-^^(rT^-^^^)=o- (4-18) 

Note that the dispersion relation (4.8) can also be cast into the form 

which facilitates a comparison with the previous equation. Expansion of Eqs. (4.18) 
and (4.19) to lowest order in Kj = kvj yields the MHD dispersion relation (4.10), 
i.e. in this case, the conservation equation (4.18) has a zero integration constant 
and simply expresses equipartition between kinetic and magnetic energy densities. 
This is not true any more if finite gyrokinetic effects are considered. 

4-2. Linear electrostatic waves 

Let us now discuss pure linear electrostatic waves, which are obtained by taking 
the trivial solutions, Vjt ~ Bt = 0, of the wave equation (3.16), which also implies 
that Ejji = 0. Then the linearized electrostatic wave equation (3.17) simply reads 
Ve En = 0. After FT, we obtain that (fc^ + X~'^)En{k) = 0, which explicitly yields 
the dispersion relation 

We may only consider here the case of zero drifts, i.e. Ujn = 0, and a simple electron- 
proton plasma. Then we always find two solutions for from the equation 

In the long- wavelength limit {k 0), the diverging phase speed VL{k) = u)p/k 
corresponds to the Langmuir wave, with the total plasma frequency being defined 
by uip = + Wp. For the ion-acoustic or sound wave, we obtain the constant speed 
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defined as 



K-c^+^gcg^ M7eTe + 7pTp)_ (4.22) 



Here fee is Boltzmann's constant. The general solution in terms of frequency follows 
from the biquadratic equation 

iJ^ - UJ^ (W^ ^ ^2 ^ (^^^)2 ^ (^^^)2^ ^ (Ce/c)2(Cpfc)2 + ^^(Cpfc)^ + a;2(c^fc)2 = Q. 

(4.23) 

In the short- wavelength limit (A; — ?■ cxd), we obtain two solutions corresponding to 
the proton-acoustic wave with w ~ kcp, or electron-acoustic wave with oj w kCg. 
Both modes are usually strongly Landau damped if a thermal Vlasov description 
of the plasma is used. The sound wave, however, can exist since We rrip, and 
thus Vs ~ i/fcB7e^e/TOp for Te > Tp, so that strong proton Landau damping can 
be avoided. 

For a multicomponent plasma with drifts, the structure of the eigenmodes be- 
comes correspondingly richer, as each species contributes its own plasma frequency 
and thermal speed, as well as drift speed. In the presence of a compressive transverse 
wave, these modes all become coupled and are driven by the nonlinear ponderomo- 
tive electric fields Ejn according to Eq. (3.17). Similarly, the transverse eigenmodes 
defined by Eq. (4.8) arc driven according to Eq. (3.16) by the longitudinal electric 
field i^n and the combined action of the various Ejn- 

4-3. Compressive Alfven-cyclotron-acoustic waves 

In this section, we will consider the coupling between the nonlinear electromagnetic 
Alfven-cyclotron waves and the electrostatic modes. We recall that no assumptions, 
such as incompressibility, had to be made as to derive the wave equations (3.16) 
and (3.17). We shall also rewrite the transverse momentum equation (3.8), which 
describes the gyromotion as a second-order wave equation. Since the longitudinal 
gyration scale Vj depends on the density according to Eq. (2.34), its differentiation 
has to be considered. If we neglect the magnetic field for a moment, then we get 
for the transverse motion an equation in the form 

dlnn," d \ „ / . „ ,^ 

-^d?j^^*=0- ^'-''^ 

Mathematically speaking, this is the well known equation for a harmonic oscillator, 
with an amplitude that may vary exponentially in ^ at a scale set by the density 
gradient length. For the differential operator yielding harmonic oscillations (first 
two terms of Eq. (4.24)), we introduce the symbol Vv^ to be used below. For the 
density gradient term, we may approximately write 




d , d , ^ , 1 d ^ Sn-j d In (5n, 



with the average constant background density fij. If the gradient is positive (neg- 
ative), and therefore the density increases (decreases), the longitudinal scale and 
thus the amplitude of Vjt will decrease (increase) correspondingly. However, as 
long as the relative density variation remains small, a few percent say, this change 
will occur on a much larger scale than rj, namely Lj = rjUj/Snj. If the density 
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fluctuates about zero, the net effect of the density modulation on Vjt will remain 
comparatively small. 

For the sake of consistency, we will, in the remainder of this section, fully retain 
this density- induced possible amplitude variation of Vjt according to Eq. (4.24), 
but later on neglect the density variations and consider all density-dependent pa- 
rameters to be fixed at their background values without denoting them explicitly 
by a barred symbol. Yet, remember that the essential and lowest-order variations of 
the densities of all species have been considered and already taken care of in Gaufi' 
law and the dynamics of E^, which indeed is of order unity as the background 
electric field is zero, and similarly in Ampere's law through the appearance of the 
electric fields and Ej^. 

We introduce, by using conserved or constant quantities, normalized variables, 
such that Vjt = Vjt/V and bt = Bt/B^. Similarly, we introduce normalized electric 
fields as follows: 

_ cEjn _ cEn 



■"' BnV BnV ^^'^^^ 

We recall Eq. (3.16), according to which Cjn is just an abbreviation for the gradient 
of a potential given by the transverse kinetic energy, for which we have in normalized 
form 

Using the same normalization for the transverse electric field (2.23), we obtain 
simply that 

et - -{h X bt), (4.28) 

which is fully determined by the solution for bt. Like et, the charge densities cTj are 
now merely auxiliary quantities and obtained from an integration of the previous 
equation (3.3), which in the new variables can be written as 

IntTj = (e„ + ejn) , (4-29) 

and then be formally integrated with the result 



aj{^) =qjnjexp 



As an outcome of all the above considerations, we can now summarize the re- 
sulting set of fiuid wave equations. Firstly, we obtain for each species' transverse 
motion a forced and amplitude-modulated harmonic oscillator equation reading 

^v, = ( ^pv^(e" + e^-)^ + -^t + X ''t 1 y • (4.31) 

Secondly, one can rewrite the mutually driven and coupled wave equations for the 
longitudinal electric field and transverse magnetic field in the concise forms 

2^sen = ^^e,n, (4.32) 
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J2 ( ^^Jt + ^2 _ y2 (en + ejn)(n X Vjt) 



(4.33) 



We recall again that up to this point of our algebraic derivations, no linearization 
has been made, and the density variations were entirely accounted for. Eqs. (4.29) 
or (4.30) permit to calculate the density of species j completely through the line 
integral over the electric fields that appear in the exponential Boltzmann factor in 
Eq. (4.30). So the density is a functional of the electric potentials. However, this 
dependence of uj on ^ may now without loss of essential physics be neglected, if the 
density fluctuations can be assumed to remain small (i.e. we do not want to con- 
sider MHD shocks (Goossens, 2003), or electrostatic shocks and double layers here). 
Thus, all scales and parameters such as Vjn,Cj,ojj, Xj,£j, and Vj, which have non- 
vanishing mean values, will, from here on, be calculated by use of the background 
number density fij, as well as the conditions for quasi-neutrality and zero longi- 
tudinal current and the center-of-momentum condition. All compressive effects are 
described by the longitudinal electric field in this approximation. Consequently, 
the nonlinear equations (4.31), (4.32), and (4.33) form a closed set, which yet will 
generally require a numerical treatment to obtain solutions. This is a task beyond 
the scope of this analytical paper. 

Note that in each of these equations, the spatial variations arc determined by the 
natural scales of the dynamics of the involved field variables, i.e. by the longitudinal 
scale rj for the transverse motions of the particles, their Debye lengths Xj for 
the charge fluctuations, respectively, their skin depths ij for the magnetic field 
penetration into the plasma driven by the transverse currents. Differential motion 
of the species j might be important and is therefore included in its drift speed 
Vjn- Its effect on the parametric instabilities has been studied in work addressing 
the modulational and decay instability of Alfven waves by considering streaming of 
alpha particles in the solar wind (HoUweg et al., 1993). If there are no differential 
motions along the mean field in the background plasma, i.e. if for all j we have 
Ujn = 0, then Vjn = —V, and thus the longitudinal gyration scale simply becomes 
rj = —V/flj, which by its definition is not for each species a positive definite 
quantity as the gyrofrequency carries the sign of the charge of the species considered. 
The factor in front of the electric field term in the above Eqs. (4.31) and (4.33) 
thus changes in the drift-free case to 1/(1 — {cj/V)"^), which becomes 1/(1 — Pj) 
for V = Va, with the species plasma beta being defined as Pj = {cj/VxY- For the 
parallel propagation considered here, only this factor contains the thermal speed, 
and therefore this factor simply becomes unity for a cold multi-species plasma 
without drifts. 



4.4- Electric field fluctuations driven by an elliptically polarized Alfven wave 

In this section, we consider the electrostatic waves, which can be generated by an 
elliptically polarized Alfven wave, which for the sake of simplicity we assume to 
be given. We shall then study the effect this wave has in generating compressive 
fluctuations driven by the spatial variation of the kinetic energy of the particles 
moving coherently in the same wave magnetic field. The starting point is Eq. (4.32), 
which can with the help of Eq. (4.27) be written as a driven oscillator equation for 
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the longitudinal electric field: 

Before we write down the wave fields, let us define a proper coordinate system. For 
the right-handed orthogonal system, we choose the unit vectors n = 63 = ei x 62, 
ei = 62 X n, and 62 = n x ei. The wave may have a wave vector k, and its 
normalized (dimensionless) magnetic field reads 

bt = biCi cos (fc^) + 6262 sin (fc^). (4.35) 

Similarly, the related flow velocity of species j is given by 

Vjt = VjiEi cos (fc^) + Vj2e2 sin (fc^). (4.36) 

The associated electric field according to Eq. (4.34) reads 

kV 

ejn = COS {kO sin (fc^ {vji - v%) . (4.37) 



On the other hand, using the original definition (3.9) of this field, we obtain the 
result 

Cjn = COS {ki) sin (A;^)(&2f ji — biVj2), (4.38) 
which by comparison of the last two equations, determines the particle velocity 
components as Vji^2 = ^2,1 ^j/(^^)- Finally, we have 

ejn = ^sm{2kO{bl-bl), (4.39) 
which can be inserted in Eq. (4.32). We obtain the forced oscillator equation 

d^en 2 sin(2fcf) 

Wc recall that the Dcbyc wave length was defined in Eq. (3.5). We used as an 
abbreviation the complex wave vector q = i/Xu, and further introduced the forcing 
amplitude as 

^ = E^y(^?-^2), (4.41) 

j ^ 

which vanishes for a circularly polarized wave with 61 — b2, hut is nonzero otherwise. 
The solution of Eq. (4.40) can be obtained by using the Green's function method 
and Fourier transformation. The convolution integral of the forcing term with the 
Green's function then yields the solution depending upon ^ in the form of another 
convolution integral, which can be calculated analytically with the (in q symmetric) 
result 

£ / sin (2fcg) + sin (gg) sin (2fcO - sin (gg) 
^''^^> ~ Wq [ 2kT~q 2k^q J ' ^^'^^^ 

Eq. (4.42) solves the original Eq. (4.40) that can easily be shown by straight- 
forward differentiation. Wc recall that by definition the square of the wave vector 
q = q{V) is given by the right-hand side of the electrostatic dispersion relation 
(4.20). Therefore, g is a real number for an appropriate choice of V. The solution is 
then related naturally with the electrostatic cigcnmodcs, i.e. the sound, ion acous- 
tic, and Langmuir waves, which we already discussed in a previous section. The 
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overall solution (4.42) apparently describes forced compressive (charge) waves, oc- 
curring as electrostatic eigenmodes, and a superposed electric wave at the second 
harmonic of the transverse Alfven pump wave, the anisotropy (due to its elliptic 
polarization) of which determines the amplitude of these driven longitudinal elec- 
tric field oscillations. L'Hopital's rule allows one to determine the behavior in the 
resonant cases. For the resonances {q ±2fc), we find 



i cos(gC) 



sin(gO 



(4.43) 



which corresponds to the amplitude of the compressible oscillation growing or de- 
caying with ^, i.e. an instability in space. 



5. Discussion and conclusions 

Starting from the multi-fluid equations of a warm plasma, we have derived and 
investigated the coupled wave equations for the particles' gyromotions about the 
mean field and for the transverse magnetic field and longitudinal electric field. It 
is a natural outcome of the electromotive forces arising from compressible Alfven- 
cyclotron waves and can be derived from a potential that is just the kinetic energy 
associated with the gyromotion in the electromagnetic wave. Electric waves are 
thus excited, which can react back on the pump wave by nonlinear effects through 
terms in its own wave equation that contains the electric field explicitly. Known 
limiting cases arc reproduced, such as the standard linear electric waves like the 
ion- acoustic or Langmuir waves of course, and for the transverse magnetic field the 
usual two branches of Alfven-cyclotron and magnetosonic-whistler waves in case 
of a two-component electron-proton plasma, or many similar related branches in 
the case of a multi-ion plasma. The main result of this paper is the closed set 
of second-order wave equations (4.31), (4.32), and (4.33), from solutions of which 
the transverse electric field and charge densities of each species can be derived 
as auxiliary quantities. To study these wave equations in more detail and to find 
their nonlinear solutions is left as a future task, which will presumably require a 
numerical treatment. 

The structure of our equations already permits to derive some qualitative conclu- 
sions and to treat some simple applications (like the effect of elliptical polarization) 
analytically. Further study is certainly required to corroborate them quantitatively. 
Apparently, the weakly compressible large-amplitude Alfven-cyclotron waves can 
drive electric fluctuations, essentially of the ion-acoustic type, along the mean field, 
and thus will naturally produce an electric field that can accelerate particles and 
will lead to heating via Landau damping in a kinetic Vlasov description. By exci- 
tation of acoustic waves, the amplitude of the driver wave will be diminished until 
a dynamic wave wave equilibrium is reached. Similar processes are clearly found in 
the direct numerical simulations (Araneda et al., 2008, 2009; Valentini et al., 2008; 
Valentini and Veltri, 2009). The third-order coupling terms in Eqs. (4.31) and (4.33) 
correspond to such three- wave processes in Fourier space, and therefore will lead to 
cascading of spectral energy and broadening of the original spectrum of the pump 
wave, which need not be monochromatic. This way a new path towards micro- and 
macro-turbulence could be opened, and a non-MHD cascade is rendered possible 
by these compressive Alfven-cyclotron-acoustic wave interactions. 
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